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ABSTRACT 



Aims. We studied the abundance of HCN, H'^CN, and HN'^C in a sample of prestellar cores, in order to search for species associated 
with high density gas. 

Methods. We used the IRAM 30m radiotelescope to observe along the major and the minor axes of L1498, L1521E, and TMC 2, 
three cores chosen on the basis of their CO depletion properties. We mapped the 7 = I — > transition of HCN, H'^CN, and HN'^C 
towards the source sample plus the 7 = 1 — > transition of N2H+ and the 7 = 2 — > 1 transition of C"*0 in TMC 2. We used 
two different radiative transfer codes, making use of recent collisional rate calculations, in order to determine more accurately the 
excitation temperature, leading to a more exact evaluation of the column densities and abundances. 

Results. We find that the optical depths of both H''CN(l-0) and HN'^^C(l-O) are non-negligible, allowing us to estimate excitation 
temperatures for these transitions in many positions in the three sources. The observed excitation temperatures are consistent with 
recent computations of the collisional rates for these species and they correlate with hydrogen column density inferred from dust 
emission. We conclude that HCN and HNC are relatively abundant in the high density zone, n(H2) ~ 10^ cm"', where CO is depleted. 
The relative abundance [HNC]/[HCN] differs from unity by at most 30% consistent with chemical expectations. The three hyperfine 
satellites of HCN(l-O) are optically thick in the regions mapped, but the profiles become increasingly skewed to the blue (LI 498 and 
TMC 2) or red (L1521E) with increasing optical depth suggesting absorption by foreground layers. 
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01. Introduction 
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Tracing the dense central regions of prestellar cores is essential 
for an understanding of their kinematics and density distribution. 
Many of the well studied nearby cores are thought to be static 
with a "Bonnor-Ebert" radial density distribution (see Bergin 
& Tafalla 2007 and references therein) comprising a roughly 
1/r^ fall off surrounding a central region of nearly constant den- 
sity. However, our understanding of the central regions with H2 
number densities of 10^ cm""* or more is limited both by our 
uncertainties about dust emissivity and molecular abundances. 
Dust emission, because it is optically thin, is a useful tracer of 
mass distribution, but both temperature gradients and gradients 
in the dust characteristics (refractive index and size distribution) 
mean that dust emission maps should be interpreted with cau- 
tion. Molecular lines offer the advantage that they permit an un- 
derstanding of the kinematics, but depletion and other chemical 
effects render them often untrustworthy. Depletion is more rapid 
at high density and low temperature and these are just the effects 
that become important in the central regions of cores. 

One interesting result from studies to date (e.g. Tafalla et 
al. 2006) is that the only molecular species to have a spatial 



* Based on observations carried out with the IRAM 30m Telescope. 
IRAM is supported by INSU/CNRS (France), MPG (Germany) and 
IGN (Spain). 



distribution similar to that of the optically thin dust emission 
are NiH^ and NH3, consisting solely of nitrogen and hydrogen. 
This is in general attributed to the claim that their abundance 
is closely linked to that of molecular nitrogen (one of the main 
gas phase repositories of nitrogen) which is highly volatile and 
hence one of the last species to condense out (Bisschop et al. 
2006). However N2 is only marginally more volatile than CO 
which is found observationally to condense out at densities of a 
few 10"* cm-^^ (Tafalla et al. 2002, 2004a) and it has thus been 
puzzling to have a scenario where C-containing species were 
completely frozen out but N2 not. This has led to a number of 
models aimed at retaining gas phase nitrogen for at least a short 
time after CO has disappeared, e.g. Flower et al. 2006, Akyilmaz 
et al. 2007, and Hily-Blant et al. 2010 (hereafter paper I). There 
have also been attempts to obtain convincing evidence for the 
existence in the gas phase of species containing both C and N at 
densities above a few 10^ to 10^ cm"-' at which CO depletes (e.g. 
Hily-Blant et al. 2008 and paper I). A partial success in this re- 
gard was obtained by Hily-Blant et al. 2008 who found that the 
intensity of '^CN(l-O) behaved in similar fashion to the dust 
emission in two of the densest cores: L183 and L1544. The im- 
plication of this is that at least some form of carbon, probably in 
the form of CO, remains in the gas phase in these objects at den- 
sities above 10^ cm"-'. Thus, a tentative conclusion is that the CO 
abundance at densities above lO*" cm"^ is much lower than the 
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canonical value of 10^^ relative to H2 (Pontoppidan 2006), but 
nevertheless sufficiently high to supply carbon for minor species. 

Putting this interpretation on a more solid foundation re- 
quires observation of possible tracers of the depleted region in 
cores carried out in a manner which will permit distinguish- 
ing species associated with high density gas from those present 
in the surrounding lower density envelope. This is complicated 
though perhaps facilitated by the gradients in molecular abun- 
dance due to depletion. While molecular abundances in general 
are expected to drop at high densities when depletion takes over, 
this is not necessarily true for all species at least within a limited 
density range (see the results for deuterated species discussed 
by Flower et al. 2006). However, proving that one is observing 
emission from high density gas involves either showing that one 
can detect transitions which require high densities to excite or 
showing that the emission comes from a compact region coinci- 
dent with the dust emission peak. The former is rendered more 
difficult by the temperature gradients believed to exist in some 
cores (see e.g. Crapsi et al. 2005). The latter requires highly sen- 
sitive high angular resolution observations. A combination of the 
two is likely to be the best strategy. 

One approach to these problems which has not been fully ex- 
ploited to date is to make use of the hyperfine splitting present 
in essentially all low lying transitions of N-containing species. 
It is clear in the first place that the relative populations of hy- 
perfine split levels of species such as HCN and CN are out of 
LTE in many situations (see Walmsley et al. 1982, Monteiro & 
Stutzki 1986, and paper I). Interpreting such anomalies requires 
accurate collisional rates between individual hyperfine split lev- 
els, but such rates can be determined from the rates between ro- 
tational levels (e.g. Monteiro and Stutzki 1986) and clearly in 
principle, they allow limits to be placed on the density of the 
emitting region. It is also the case that species of relatively low 
abundance like H'^CN are found to have rather minor deviations 
from LTE between hyperfine levels of the same rotational level 
(see paper I) and in this case, one can determine the line excita- 
tion temperature and optical depth based on the relative hyper- 
fine satellite intensities. While this is clearly doubtful procedure 
in that the non-LTE anomalies can cause errors in the inferred 
optical depth, it is nevertheless (as we shall discuss) an approach 
which can give a zero order estimate of relative rotational level 
populations and hence of the local density. 

In the case of more abundant species such as the main iso- 
topologues of HCN or CN, one finds that "self absorption" or ab- 
sorption in foreground relatively low density material can oblit- 
erate any signal from the high density central regions of a core. 
However, from these tracers, one can in principle glean informa- 
tion on the kinematics of foreground layers of the core. Thus, 
one can hope to find evidence for infall. We indeed find in this 
study that profiles change in gradual fashion as a function of 
transition line strength and this does indeed yield evidence for 
either infall or expansion of foreground layers. 

Here we extend the results of paper I to three other cores 
which have been chosen on the basis of their CO depletion 
properties (see e.g. Bacmann et al. 2002, Brady Ford & Shirley 
2011). We in particular chose objects with large CO depletion 
holes indicative of a relatively large age. This might occur for 
example in sources where magnetic field is capable of slowing 
down or preventing collapse. We also used as a guide the an- 
gular size measured in the N2H^(l-0) transition by Caselli et 
al. (2002b) which is thought to correspond roughly to the area 
of CO depletion. Thus, we chose the cores L1498 and TMC 2 
which appear to have large CO depletion holes (0.08 parsec in 
the case of TMC 2). As a comparison source, we included in our 



study LI 52 IE, which shows little or no depletion in C-bearing 
species (Tafalla & Santiago 2004b). 

In this paper then, we present IRAM 30m observations of 
the / = 1 ^ transition of HCN, H'^CN and HN'^^C along the 
major and the minor axes of the three selected sources, L1498, 
L1521EandTMC2,as wellasthey = 1 ^ transition of N2H+ 
and the y = 2 ^ 1 transition of C'^O towards TMC 2. In Sect. 
2, we discuss the observational and data reduction procedures, 
summarising the results in Sect. 3. In Sect. 4 we compare the 
distribution of the line emission with respect to the dust emis- 
sion, investigating about the presence of depletion. In Sect. 5 
we describe the evidence of the asymmetry of the HCN(l-O) 
line profiles. In Sect. 6 we explain the methods used for de- 
termining the excitation temperature, while in Sect. 7 we give 
the estimated column densities and abundances for the observed 
lines. In Sect. 8 we give our conclusions. Comment on non- 
LTE hyperfine populations are provided in Appendix A, while 
in Appendix B we summarise the spectroscopic data and obser- 
vational parameters of the observed molecules. 

2. Observations 

L1498, L1521E, and TMC 2 were observed in July 2008, us- 
ing frequency switching and raster mode, with a spacing of 20", 
25" and 50", respectively, along the major and the minor axes. 
The axes were identified from the continuum emission maps 
showed in Tafalla et al. (2002) for L1498, Tafalla et al. (2004a) 
for L1521E, and Crapsi et al. (2005) for TMC 2. Observing pa- 
rameters are given in Table B.3. 

Observations of the 89 GHz HCN(l-O) and the 86 GHz 
H'^^CN(l-O) multiplets (see Table B.l) plus the 87 GHz 
HN'''C(1-0) multiplet (see Table B.2) were carried out using 
the VESPA autocorrelator with 20 kHz channel spacing (cor- 
responding to about 0.069 km s"') with 40 MHz bandwidth 
for HCN(l-O) and 20 MHz bandwidth for Hi3CN(l-0) and 
HN'^C(l-O). The final rms, in T^^, unit is crj ~ 50 mK, 
~ 20 mK, and ~ 30 mK for HCN(l-O), H'^^CN(l-O), and 
HN'^C(l-O), respectively. 

Observations of the 93 GHz N2H+(l-0) multiplet and of the 
219 GHz C'**0(2-l) Une were carried out in TMC 2 using 10 
kHz channel spacing with 40 MHz bandwidth for N2H+(l-0) 
and 20 kHz channel spacing with 40 MHz bandwidth for 
C^^O(2-l). For N2H^ the final rms, in T^\, units, in channels of 
width 5v = 0.031 km s-i, is o-T ~ 100 mK while for C^'*0(2-l), 
the final rms (channels of width 5v = 0.027 km s"') is o-j ~ 250 
mK. 

Data reduction and analysis were completed using the 
CLASS software of the GILD AS' facility, developed at the 
IRAM and the Observatoire de Grenoble. In what follows, all 
temperatures are on the main-beam scale, T,^b - Ff.ffT*^/Bfiff, 
where is the antenna temperature corrected for atmospheric 
absorption, while Beff and Feff are the beam and the forward effi- 
ciencies, respectively (see Table B.3 for the numerical values of 
efficiencies). 

3. Observational results 

Figure 1 shows the continuum emission of the three cores to- 
gether with the positions mapped in lines described in Sect. 2. 
The selected sources have already been widely observed in the 
past as possible tracers of conditions in the high density core 

' http://www.iram.fr/IRAMFR/GILDAS 
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L1498 L1521E TMC 2 




Aa [arcsec] Aa [arcsec] Aa [arcsec] 

Fig. 1. Observed positions in L1498 (left panel), L1521E (middle panel), and TMC 2 (right panel), superposed on the dust emission map smoothed 
to 28", from Tafalla et al. (2002), Tafalla et al. (2004a), and Crapsi et al. (2005), respectively. Mapped positions in HCN(l-O) {white crosses), 
H'''CN(l-0) and HN'3C(l-0) (white squares), and N2H+(l-0) (yellow circles). Contours represent 35, 50, 65, 80 and 95 per cent of the peak 
value of the dust emission which is 18.0, 21.0, and 45.0 mjy/(ll" beam) for L1498, L1521E, and TMC 2, respectively. The (0,0) position 
corresponds to q'(2000) = 04'>10™51.5% (5(2000) = 25°09'58" for L1498, to a(2000) = 04''29™15.7% 5(2000) = 26°14'05" for L1521E, and to 
^(2000) = 04*'32"48.7% <5(2000) = 24°25'12" for TMC 2. 



nucleus. We succeeded in observing HN^C(l-O) and the three 
hyperfine components of HCN(l-O) and Hi3CN(l-0). 

Figure 2 shows a comparison of the line profiles for the dif- 
ferent tracers. We plotted the weakest (F = — > 1) component 
of HCN(l-O) at 88633.936 MHz, the strongest (F = 2 ^ 1) 
component of H'3CN(l-0) at 86340.184 MHz, HN'^Ql-O) at 
87090.675 MHz, the isolated component (Fi,F = 0, 1 -> 1,2) 
of N2H+(l-0) at 93176.2650 MHz, and the C'^0(2-l) transi- 
tion at 219560.319 MHz. We used the components enumerated 
here also for the comparison between line and dust emission de- 
scribed below, supposing them to be the most optically thin. The 
offsets for each source are relative to the dust peak emission (see 
caption of Fig. 1). 

There are clear trends in the line widths which we derive for 
our three sources with values of order 0.2 km s"' for L1498, 0.3 
km s"' for L1521E, and 0.45 km s"' for TMC 2. Besides, we no- 
ticed that there is a reasonable agreement between the systemic 
velocity, Vlsr, of N2H+ and H'^CN to within 0.08 km s while 
line widths of H'^'CN seem larger by 0.05 km s"' with respect to 
N2H^, but this may be due to '''C hyperfine splitting (Schmid- 
Burgk et al. 2004). 

There is a clear difference between the HCN and H^^CN pro- 
files in the three sources: while in L1521E, the HCN profile is 
broader than its isotopologue, in L1498 and TMC 2, HCN and 
H^^CN have very different profiles, presumably because of high 
optical depth in HCN. 

All these molecules show hyperfine structure, although in the 
case of HN'^C it is not possible to avoid the blending of the 
components, and the hyperfine fitting has been conducted using 
the HFS method in CLASS. To fit this line, we followed van 
der Tak et al. (2009) who found that it consists of eleven hyper- 
fine components which can however be reduced to four "effec- 
tive" components. As an instance, in the upper panel of Fig. 3, 
we show the fit of the HN'^C Hne at the offset (10, 20) towards 
L1498 together with the four distinguishable hyperfine compo- 
nents. The fit gives a value for the total optical depth equal to 
T = 4.58 + 0.32 and all the observed points show similar val- 



L1498 (10",20") L1521E (-25",0") TMC 2 (-50",0") 




Vlsr [km s"'] V^sh [km s '] Vi^r [km s"'] 

Fig. 2. Emission line for the different tracers of the source sam- 
ple towards the dust emission peak. The 88633 MHz component of 
HCN(l-O) is shown together with the 86340 MHz component of 
H"CN(l-0), HN'3C(l-0), the 93176 MHz isolated component of 
N2H+(l-0), and C'*0(2-l). Spectra have been multiplied by a scal- 
ing factor to allow the simultaneous comparison. Red dotted lines show 
the systemic LSR velocity of the sources evaluated from the well- 
determined frequencies of N2H+(l-0), see Pagani et al. (2009). 



ues, as shown in Figure 4, with an average of (r) = 5.16 + 0.86. 
Also for H'^CN we found a good simultaneous fit of the three 
hyperfine components (see lower panel of Fig. 3) and, as HN'^C, 
this line is somewhat optically thick with a total optical depth of 
T - 3.11 + 0.77 at the offset (10,20) and a mean optical depth 
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Fig. 3. HN"C(l-0) (upper panel) and H'3CN(l-0) (lower panel) emis- 
sion towards L1498 at the offset (10,20); black solid histogram, hy- 
perfine components used for the fit, blue dashed lines, centred at the 
frequencies listed in Table B.2, and line fit, magenta solid line. 
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Fig.4. Total optical depth of HN'^C(l-O), solid squares, and 
H"CN(l-0), empty squares, for all the observed positions in L1498 
as a function of the molecular hydrogen column density. 



of (t) = 4.34 + 1.15. Hence one can expect the optical depth in 
HCN hyperfine components to be of order 30-100 assuming the 
canonical isotopic ratio for ['2C]/['^C] of 68 (Milam et al. 2005). 



4. Probing the presence of depletion 

In this section, we study the dependence of the integrated in- 
tensity, W, of the HCN, H'^'CN and HN'^'C transitions on offset 



from the dust emission peaks. In the case of HCN, we consider 
the weakest F = — > 1 line and in the case of H'^CN, the 
strongest F = 2 — > 1 satellite. The observational results are 
shown in Figures 6, 7, and 8 where we compare the observed 
intensity of these lines with the dust emission which is taken to 
be representative of the Ht column density. It is useful first to 
comment briefly on what one can expect qualitatively to learn 
from such studies. 

One notes in the first place that optically thick transitions, 
due to scattering and line saturation, are naturally likely to have 
a broader spatial distribution than thin (or very moderate optical 
depth) transitions. Thus, the HCN(F = — » 1) transition can 
be expected to be roughly an order of magnitude more optically 
thick than H'^CN(F = 2^1) assuming the local interstellar 
['^C]/['^C] ratio (Milam et al. 2005) and no fractionation. One 
thus expects a broader spatial distribution of HCN than H'-'CN 
and indeed this is confirmed by the results shown in all the three 
sources. For example, the half power size of HCN in L1498 
along the NW-SE cut is 170" as compared with 120" for H'^CN 
and HN'-'C, and 180" for the dust emission. One concludes that 
the '^'C substituted isotopologue is a much better tracer of HCN 
column density than the more abundant form even when using 
the weakest hyperfine component. 

However, one also notes from these figures that the C'^O dis- 
tribution is essentially flat as has been seen in a variety of studies 
(Tafalla et al. 2002, 2004a). This has been attributed to depletion 
of CO at densities «(H2) above a critical value of order a few 
times 10"* cm""* and our observations are entirely consistent with 
this. It is also true however (see also Tafalla et al. 2006 and pa- 
per I) that species like H'^CN and HN'^'C, while they clearly 
have different spatial distributions from the dust emission, have 
half power sizes which are roughly similar (see results for L1498 
above). Such effects clearly do not apply to NaH^ which does 
not contain carbon and which in general has been found to have 
a spatial distribution similar to that of the dust (Tafalla et al. 
2006). However, one sees in Fig. 6-8 that HN'^^C and H'^CN 
have a spatial distribution closer to that of the dust emission and 
we conclude that these two isotopologues may often be very use- 
ful tracers of kinematics as N2H+ at densities above the critical 
values at which CO depletes. It is possible that CO, while de- 
pleted by roughly an order of magnitude relative to the canonical 
[CO]/[H2] abundance ratio of lO""*, is nevertheless sufficiently 
abundant to account for minor species such as HCN and HNC. 

Finally, we note from Fig. 8 that towards TMC 2, where 
the dust peak H2 column density is larger than in the other two 
sources, there seems to be a good general accord between the 
dust emission and the intensity of HN'^'C, H'^CN, and N2H+. 
The CO however has again a flat distribution suggesting once 
more that it is depleted in the vicinity of the dust peak. There 
is no reason in this source to suppose that N2H^ traces the high 
density gas around the core nucleus better than the '^C substi- 
tuted isotopologues of HCN and HNC. There are however slight 
differences close to the dust peak between the different species 
which are perhaps attributable to excitation and optical depth 
effects, but could also be caused by the depletion which presum- 
ably all molecules undergo if the density is sufficiently high. We 
conclude that cyanides and isocyanides are useful tracers of gas 
with densities around 10^ cm"-'. 

We note that models of a collapsing prestellar core indeed 
suggest that HCN and HNC should remain with appreciable gas- 
phase abundance at an epoch when the CO abundance has de- 
pleted to a few percent of the canonical value of 10"^ relative 
to H2. To illustrate the differential freeze out of HCN relative to 
that of CO found in the chemical models. Fig. 5 displays the de- 
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Fig. 5. Tlie HCN depletion as a function of the CO depletion during the 
collapse, computed by the gravitational collapse model of Hily-Blant 
et al. (2010), see their Fig. 8. Values are normalised with respect to the 
initial abundances. Each point is labelled with the time of the collapsing 
gas. Two models are shown: a free fall model labelled % (red dashed 
cun'e) and a collapsing model with a time scale multiplied by a factor 
ten, labelled 10f(f (blue solid curve). The black dotted cun'e shows the 
positions of equal depletion for CO and HCN. 



pletion of HCN as a function of the depletion of CO computed 
in the collapsing gas, starting from steady-state abundances at a 
density of 10"* cm"^. The numerical gravitational collapse model 
is described in Sect. 6.2 of Hily-Blant et al. (2010). Two mod- 
els are shown for two different collapse time scale: free-fall time 
(fff) and lOfff. It may be seen from Fig. 5 that the depletion of 
HCN relative to that of CO strongly depends on the collapse 
time scale; the model shows that HCN is depleted on longer 
timescales than CO, in particular if the dynamical timescale is 
several free fall times. HNC, not shown in the plot, has the same 
behaviour of HCN. 



5. Line profiles 

In this Section, we discuss the behaviour of the HCN line profiles 
as a function of optical depth. It is clear from Fig. 2 that towards 
L1498 and TMC 2, the F = ^ 1 component of HCN has 
a profile skewed to the blue relative to F = 2 — > 1 of H'^CN 
and it is tempting to interpret this as being due to absorption in 
a foreground infalling layer. Fig. 9 illustrates this well showing 
all three HCN components compared with the strongest F - 
2 — > 1 component of H'^'CN. One sees that there is in fact a 
progression going from the highly skewed optically thick F - 
2 —> 1 component of HCN to the essentially symmetric profile 
of H'^CN. In the following, we attempt to quantify this trend. 

Looking at Fig. 9, related to the three hyperfine compo- 
nents of HCN and the H'''CN(F - 2 —> I) component, one 
can qualitatively argue that the greater the relative intensity of a 
line, the higher is the skewness degree. It can be seen that the 
H'-^CN(F = 2^1) component (gray histogram) is fairly sym- 
metric and indeed, as seen earlier, its optical depth is not high 
since it is next to optically thin limit, but HCN components are 
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Fig. 6. Comparison between the dust emission (gray histograms) and 
the integrated intensity of the observed species along the two cuts 
in L1498. Upper panels: solid squares, HCN(l-O) [x3]; empty cir- 
cles, C'**O(l-0) [xO.8]; solid circles, C"*0(2-l) [xl.4]. Lower pan- 
els: empty triangles, H'^CN(l-O) [xlO]; empty squares, HN"C(l-0) 
[x2.5]; empty pentagons, N2H^(l-0) [x3]. The typical error on the in- 
tegrated intensity is about 20 mK km s"' . 
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Fig. 7. Comparison between the dust emission (gray histograms) and 
the integrated intensity of the observed species along the two cuts in 
L1521E. Upper panels: solid squares, HCN(l-O) [x2.5]; solid circles, 
C"*0(2-l) [xO.4]. Lower panels: empty triangles, H'''CN(l-0) [x27]; 
empty squares HN'^C(l-O) [x5.5]; empty pentagons, N2H+(l-0) 
[xlO]. The typical error on the integrated intensity is about 20 mK km 



skewed towards the blue in L1498 and TMC 2, while in L1521E 
they are skewed towards the red. 

The superposition of the different hyperfine components of 
HCN led us to the discovery of a correlation between the line 
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Fig. 8. Comparison between tlie dust emission (gray histograms) and 
the integrated intensity of the observed species along the two cuts in 
TMC 2. Upper panels: solid squares, HCN(l-O) [x4]; solid circles, 
C'*0(2-l). Lower panels: empty triangles, H''CN(l-0) [x25]; empty 
squares, HN'^CCl-O) [x3.5]; empty pentagons, N2H+(l-0) [x4.5]. The 
typical error on the integrated intensity is about 20 mK km s"' . 



profile and its intensity. In particular, a red-absorbed line profile 
(that is skewed towards the blue) is a hint for the presence of 
an outer layer which is absorbing the emission of the inner layer 
while moving away from the observer, suggesting infall motions. 
Conversely, a blue-absorbed line profile (that is skewed towards 
the red) is an indication of the motion of the outer absorbing 
layer towards the observer, that is outflow motions. 

To quantify these deviations from the expected Gaussian 
shape in the optically thin limit, we compared the asymmetry 
degree towards the different positions mapped, using the defini- 
tion of skewness, 6V, given by Mardones et al. (1997) 



6V^ 



V, 



thick 



-V,y 



AVthin 



(1) 



where Vthick and Vthm are the velocities at the peak of the op- 
tically thick and thin components, respectively, while AVthin is 
the line width of the thin component. The normalisation of the 
velocity difference with AVthin reduces bias arising from lines 
of different width measured in different sources, allowing a 
more realistic comparison of the values of dV in our sample. 
Velocities and line widths have been determined by Gaussian 
fits; in the event that optically-thick line profiles were double 
peaked, we fitted these two components with two Gaussians, as- 
signing to Vthick the central velocity of the Gaussian relative to 
the stronger of the two peaks. We supposed the F - 2 ^ I com- 
ponent of H'^'CN to be optically thin, so that Vhin = V'hi'cn ™d 
AVthin - AVhucn- a value of 6V lower than zero suggests a red 
absorption of the line, while a positive value a blue absorption. 

Fig. 10 shows the values of the skewness degree, 6V, for the 
three hyperfine components of HCN as a function of the column 
density of molecular hydrogen for our source sample. Hence we 
confirm the suggestion made earlier that: (;) the absolute value 
of 6V is greater for the strongest hyperfine components, see also 
Fig. 11, in fact, in all the three sources, \6V{F = 2 — » 1)| > 



\dV(F = 1 1)1 > \6V{F = ^ 1)1; (;7) in L1498 and TMC 2, 
the emission lines are red-absorbed, since 6V < 0, and this is a 
hint for the presence of infall motions; (;//) in L1521E, the emis- 
sion lines are blue-absorbed (6V > 0), suggesting expansion; (iv) 
as expected, 6V decreases from the centre to the outer part of the 
cores, that is for decreasing values of A^(H2), dropping with line 
intensities, though sometimes not much; (v) 6V{F = 0^1) for 
HCN seems to be rather independent of the H2 column density, 
being almost constant for L1498 and L1521E, probably because 
this hyperfine component is the weakest line and it is closer to 
the optically thin limit. 



6. Excitation temperature results 

Even though real deviations from LTE populations are present 
(see Appendix A), the LTE assumption is an approximation 
which is useful for many purposes and we calculated the ex- 
citation temperature, Tex, from a simultaneous LTE fit of the hy- 
perfine components in the observed species, using the measured 
intensity of optically thick transitions. The main-beam tempera- 
ture is given by 



(2) 



where Jv{T) - To/[exp(To/T) - 1] is the Planck-corrected 
brightness temperature and Tbb - 2.73 K is the temperature of 
the cosmic background. Tq = hv/k, where v is the transition 
frequency, and li and k represent Planck's and Boltzmann's con- 
stants, respectively. If both the source and the beam are Gaussian 
shaped, the beam filling factor, /b, is given by 



/b = 



fls 



Qb + ' 



(3) 



where Q.b - 1.1330g and Us = 1.1 330^ denote the solid angles 
covered by the beam and the source, respectively, while 6^ and 
6s are the half-power beamwidth of the beam and the source, 
respectively, this latter evaluated by taking the line intensities 
stronger than the 50% of the peak value. We found /b to be equal 
to unity for most of the tracers, because Qs is at least one order 
of magnitude greater than Qb, except for H'^CN and HN'^C in 
L1521E(/b = 0.88). 

The correct determination of the excitation temperature for 
HN'^'C and H'^CN is fundamental. In fact, the gas density of 
the source sample (~ lO*" cm"-', see Tafalla et al. 2004a, Tafalla 
& Santiago 2004b, and Crapsi et al. 2005 for L1498, L1521E, 
and TMC 2, respectively) corresponds to a region where the Tgx 
of the two isotopologues changes rapidly with density, namely 
between the radiation- and collision-dominated density ranges. 
We can thus in principle use our Tgx determinations to constrain 
the density and temperature in the region where the H'^'CN and 
HN'^C lines are formed. Since these transitions are emitted in 
the region where CO is depleted (see Sect. 4), we gain informa- 
tion about the high density "core of the core". It is worth noting 
that our observations also are a test of the collisional rates for 
HCN and HNC (Sarrasin et al. 2010, Dumouchel et al. 2010). 
We illustrate this in the next section. 



6.1. Comparison of observed and expected excitation 
temperature 

In Fig. 12, we show the comparison between the values of the 
excitation temperatures of H'^CN and HN'^C evaluated from 
the simultaneous fit of the hyperfine components for the three 
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Fig. 9. Superposition of tlie tliree hyperfine components of HCN( 1 -0) (F = 1 ^ 1 , red line; F = 2 — > 1 , green line; F = — > 1 , blue line) with the 
strongest component (F = 2 — » 1, gray histograms) of H''CN(l-0) [x4], see Table B.l for component indices. In order to compare line shapes 
and intensities, components have been shifted in frequency. 
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Fig. 10. Degree of skewness as a function of the column density of molecular hydrogen: Si = (V, - Vh'3cn)/'^^h"cn> where i = 1,2,3 is the 
ith component of HCN(l-O), and VhIj^n and AVhijcn are the velocity and the line width of component 2 (F = 2 — > 1, the strongest) of the 
isotopologue H'^^CN(l-O). Values of SV from observations of component 1 (F = 1 — > 1), triangles; component 2 (F = 2 — » 1), squares; 
component 3 (F = — > 1), pentagons. See also Table B.l for component indices. 
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Fig. 11. Degree of skewness as a function of the relative intensities of the hyperfine components of HCN(l-O). Linear regressions emphasise 
the relationship between SV and line strength {dashed red lines). Values of SV from observations of the F = 1 ^ 1 component {triangles), the 
F = 2 — > 1 component {squares), and the F = — » 1 component {pentagons). 



sources examined, revealing that rex(HN'^C) is essentially al- 
ways greater than TexCH'-'CN) which we presume to be due to 
the differing collisional rates. In the same plot, two curves trace 
the values for Tex computed using RADEX (van der Tak et al. 
2007) which uses the new collisional rates for HNC (Sarrasin 
et al. 2010, Dumouchel et al. 2010) that are no longer assumed 
equal to HCN. These curves assume kinetic temperatures, Tkin, 
equal to 6 and 10 K, and typical values of 2 x 10'" cm"- and 0.2 
km s ' for column density and line width, respectively. We also 
assumed that the species are cospatial and hence come from re- 



gions of the same kinetic temperature. Notice that we rely on the 
fact that these lines are not too optically thick and thus Tex de- 
terminations are relatively insensitive to column density and line 
width. We notice an impressive good agreement between theory 
and observations, even if L1498 shows values a little below the 
expected theoretical trend. 

For L1498, we conclude from the RADEX results that the 
observed Tex are consistent with a density of a few times 10"* 
cm"-' but this estimate is sensitive to the assumed temperature. 
If the kinetic temperature is 10 K, then we can exclude densi- 
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Fig. 12. Excitation temperatures of H"CN(l-0) versus HN'^C(l-O) 
evaluated from observations of L1498 (yellow squares), L1521E {green 
circles), and TMC 2 (black triangles). The solid and dashed black lines 
trace the values for T^x computed using RADEX for kinetic tempera- 
tures, r^n, of 6 and 10 K; red empty circles and red solid circles show 
temperatures where number density assumes values equal to 10* and 
10^ cm"^^ for = 6 K and 10 K, respectively. The dotted black line 
depicts the positions where the two excitation temperatures would be 
equal. 



ties as high as 10^ cm"^ in the region where H'^CN and HN'^C 
are undepleted. On the other hand, for a temperature of 6 K, 
such densities are possible. Clearly, more sophisticated models 
are needed to break this degeneracy. In general, our results are 
consistent with the conclusion of Tafalla et al. (2006) that the 
HCN-HNC "hole" is smaller than that seen in CO isotopologues. 
This suggests that some CO, but much less than the canonical 
abundance of about 10 is present in the CO depleted region to 
supply carbon for HCN and HNC. 

For LI 52 IE and TMC 2, a central kinetic temperature of 6 
K seems possible. From dust emission models, they are thought 
to reach central densities of around 3 x 10'' cm"^ (Tafalla & 
Santiago 2004b, Crapsi et al. 2005) which would correspond 
to excitation temperatures for H'^CN and HN'^C higher than 
observed if the temperature was as high as 10 K. 

In Fig. 13, we show the observed excitation temperature of 
HN'^C and H'^'CN against the H2 column density inferred from 
the dust emission. In all the three sources, one notes a clear cor- 
relation between Tex and A^(H2). This is an important confirma- 
tion of the hypothesis that Tex is a measure of collisional rates 
and that the dust emission peak is also a peak in hydrogen num- 
ber density. It is also noteworthy that in TMC 2, which has the 
largest central column density of our three sources, the excita- 
tion temperatures continue to increase up to column densities of 
8 X 10^2 cm 2. This suggests that HN^^C and H'^CN are still 
present in the gas phase at densities close to the maximum value 
in TMC 2 (3 X lO'^ cm"^ according to Crapsi et al. 2005). On 
the other hand, the Tex values at the dust peak of TMC 2 are not 
greatly different from that measured at the dust peak of L1498 
suggesting perhaps that the central temperature is lower in TMC 
2. In any case, we conclude that excitation temperature deter- 
minations of H'^CN and HN'^C are a powerful technique for 
investigating physical parameters towards the density peaks of 
cores. 



6.2. Monte Carlo treatment of radiative transfer in L1498 

L1498 has been already extensively studied in molecular lines 
(Tafalla et al. 2004a, 2006, Padovani et al. 2009) and we used 
the Monte Carlo radiative transfer code in Tafalla et al. (2002) to 
model H'^CN and HN'^C exploiting the recent collisional rate 
calculations (Sarrasin et al. 2010, Dumouchel et al. 2010). The 
core model is the same as the one used in the analysis of the 
molecular survey of L1498 published in Tafalla et al. (2006). The 
core has a density distribution derived from the dust continuum 
map and an assumed constant gas temperature of 10 K, as sug- 
gested by the ammonia analysis. The radiative transfer is solved 
with a slightly modified version of the Monte Carlo model from 
Bernes (1979). The molecular parameters were taken from the 
LAMDA- database, where the rates are computed for collisions 
with He and we assumed that the H2 rates are larger than the He 
rates by a factor of 2, the same criterion being used for HCN in 
Tafalla et al. (2006). 

With regard to H'^CN, Tafalla et al. (2006) made use of 
the collisional rates of Monteiro & Stutzki (1986) which con- 
sider hyperfine structure complemented for higher transitions 
with those of Green & Thaddeus (1974) for higher J. In this 
model, the abundance law is a step function with an outer value 
of 1.7 X 10"'° and a central hole of 8 x 10"' cm. We ran a 
model with the same abundance law, but this time using the 
new LAMDA file for HCN, which does not account for hyper- 
fine structure. To compensate for the lack of hyperfine structure, 
which spreads the photons in velocity, we broadened the line by 
increasing the amount of turbulence from 0.075 km s"' to 0.2 
km s"' . The upper panel of Fig. 14 shows this alternative model 
which nicely fits the radial profile of observed H'^^CN intensities 
and this allows us to conclude that the use of the LAMDA file 
for HCN without hyperfine structure has little effect on the abun- 
dance determination, the result being as good as that of Tafalla 
et al. (2006). 

For HN'^'C, we used the collisional rates for the main iso- 
topologue, HNC. As for H'^'CN, the HNC hyperfine structure is 
not taken into account, so we set again the parameter relative to 
the turbulent velocity to 0.2 km s"'. The lower panel of Fig. 14 
shows the observed HN'^C integrated intensities for two fits: one 
assumes equal abundances for H'^CN and HN'^C, the other as- 
sumes that the HN'-'C abundance is 1.2 times that of H'^^CN. 
These radial profiles suggest an abundance ratio between the two 
molecules close to 1, with a likelihood of being closer to 1.2 to- 
wards the centre. 

The Monte Carlo code also computes the radial profile of 
excitation temperature, which is presented in the upper panel of 
Fig. 15. We notice that Tex for both species decreases with the ra- 
dius from about 4.8 K and 3.8 K for HN'^^C and H'^^CN, respec- 
tively, in the core interior, towards the cosmic background tem- 
perature near the outer edge of the core. Even more interesting, 
we found rex(HN'^C) systematically higher than rex(H'^CN), 
in nice agreement with the values computed from observations 
(see lower panel of Fig. 15 compared with Fig. 12). It is im- 
portant to stress that Tex values from observations refer to some 
kind of weighted mean along the line of sight for different po- 
sitions, while Tex values from Monte Carlo analysis are related 
to each shell in the core. This means that the two numbers are 
closely connected, but do not exactly have the same meaning 
and a proper comparison would require simulating the hyperfine 
analysis in the Monte Carlo spectra. However, the effect is global 
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Fig. 13. Excitation temperatures of H'^^CN(l-O) and HN''C(l-0), empty and solid points, respectively, as a function of the molecular hydrogen 
column density from observations of L1498 {yellow squares), L1521E {green circles), and TMC 2 {black triangles). 
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Fig. 14. Radial profile of observed H"CN(l-0) and HN"C(l-0) in- 
tegrated intensities (upper and lower panel, respectively) and model 
prediction for a core with a central hole of 8 x 10'* cm and an outer 
abundance value of 1.7 x 10"'" {red solid lines) and 2.04 x 10"'" {blue 
dashed line). HN'^CCl-O) data come from this study and H'^CNCl-O) 
data come from this study and Tafalla et al. (2006). 



over the cloud and affects every layer, so it cannot be ignored in 
any analysis. 



7. Column density and abundance estimates 

Based on our earlier discussion, HCN hyperfine components are 
too optically thick to be used for column density determinations. 
This means that, while H'^CN, as well as HN'^^C, originate in 
the central part of the core, HCN is dominated by the foreground 
layer emission and hence it is not possible to observe the centre 
of the core. Besides, there is an unequivocal difference between 
the Hi^CN(F = 2^1) and the HCN(F = 0^1) profiles 
(see Fig. 9). In fact, H'-'CN seems to be much more Gaussian 
than the weakest HCN line and it is likely to underestimate the 
optical depth of HCN. 





T^^(HN C) [K] 

Fig. 15. Upper panel: radial profile of excitation temperature for 
H'3CN(l-0), red solid line, and HN'''C(l-0), blue dashed line, in 
L1498 as predicted by our best fit Monte Carlo model; the vertical 
dashed line represents the central hole radius and the horizontal dashed 
line shows the cosmic background temperature limit. Lower panel: 
comparison between HN"C(l-0) and H'^^CN(l-O) local values of the 
excitation temperature in each core shell; the dashed line depicts the 
positions where the two excitation temperatures would be equal. 



However, we found that even HN'-^C and H'^CN are moder- 
ately optically thick (see Sect. 3) and in this case one can derive 
the column density in the y'th level, Nj, integrating Eq. (2) over 
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frequency to obtain the integrated intensity, Wj, with the optical 
depth given by 

T = -^AjiNjie''''-"'^- - 1 ) 0( V) , (4) 

where Aji is the Einstein coefficient, and (p{v) the profile func- 
tion, which is a sum of Gaussians (assuming a maxwellian distri- 
bution of the particle velocities) with the appropriate weights and 
shifts with respect to the central frequency, properly accounting 
for the hyperfine structure. As shown in Fig. 16, for small optical 
depths (t < 1) Nj is in direct ratio to Wj while, as t increases, 
the curve flattens even if the flattening is not very sharp because 
when the main component is thick, the satellites are still thin. 
Allowing for increasing optical depth is important not to under- 
estimate column density: as an instance, at the emission peak 
of H'^CN in L1498, where t = 4.58 (see Sect. 3), the linear 
approximation causes an inaccuracy of about 45% and this per- 
centage mainly depends on the optical depth, but not on Tex. In 
fact, for decreasing excitation temperatures, the linear approxi- 
mation decreases its slope, but also the flattening due to opacity 
is accentuated. HN'^C indicates a similar deviation. 







1 1 1 1 1 1 1 1 1 1 
linear 




1 1 1 1 1 1 1 1 1 1 _ 






approximation 52%^ 




-^^T = 6.0 




0,6 


\ 47%/- 




-^>J<S5.0 






\ 42% / - - 




>;<=4.o 






35%/ 
















B 


0.4 












26% / -i<^=2.0 
























flattening due 




0.2 






to opacity 




0.0 


- /t = 







0.0 0.5 1.0 1.5 2.0 2.5 

Nj [10'^ cm"^] 



Fig. 16. H''CN(1-0) integrated intensity as a function of the column 
density of the 7=1 level for = 3.56 K (the highest excitation tem- 
perature observed in L1498). The percentage values show the deviation 
of the linear approximation with respect to the correct value taking ac- 
count for optical depth. 



Hence, we determine the column density Nj corresponding 
to the observed W j carrying out this procedure using the max- 
imum and the minimum value of the excitation temperature in 
order to estimate errors in A^^. It is important to remark that in 
this way column densities are directly estimated from the inte- 
grated intensity of the spectra, avoiding the use of optical depth 
values evaluated from the fit of the hyperfine components which 
have often large uncertainties (see Fig. 4). Finally, for the total 
column density, A^, it holds 

N Q ' 

where Q - "^J^^Qgj e^^'l'^^" is the partition function, Ej and gj 
being the energy and the statistical weight of the upper yth level, 
respectively. 



Hence, we derive the column density and the abundance of 
these species with respect to molecular hydrogen for the main 
isotopologues HNC and HCN, assuming the canonical isotopic 
[i2C]/[i3C] ratio (Milam et al. 2005). Figure 17 shows HNC and 
HCN column densities as a function of H2 column density which 
has been calculated from the millimetre dust emission, assuming 
a dust opacity per unit mass of 0.005 cm- g"' and a dust temper- 
ature of 10 K. Notice the correlation between HNC and HCN 
column densities and their optical depths in Fig. 4 for L1498: 
for A^(H2) > 2.5 x 10^- cm^^, A?(HNC) becomes higher than 
A^(HCN) as wen as t(HNC) > t(HCN). Besides, HNC and HCN 
abundances increase towards the peak in TMC 2 but not in the 
other two cores. 

As an instance, in Table 1 we show the values of column 
density and abundance with respect to H2 for HNC and HCN at 
the dust emission peak. 

Table 1. Column densities and abundances for HNC and HCN at the 
dust emission peak in our source sample. 



source 


yvtHNC) 


A'(HCN) 




[10'-* cm--] 


[10'-* cm-2] 


L1498 


1.58(0.56) 


1.24(0.18) 


L1521E 


0.51(0.14) 


0.47(0.12) 


TMC 2 


2.41(0.42) 


1.86(0.26) 




[HNC]/[H,] 


[HCN]/[H2] 




[10-'] 


[10-'] 


L1498 


5.00(1.95) 


3.92(0.96) 


L1521E 


1.36(0.42) 


1.27(0.37) 


TMC 2 


3.05(0.57) 


2.35(0.36) 



From the point of view of the chemistry, it is even more in- 
teresting to check the ratio of abundances [HNC]/[HCN] as a 
function of the column density of molecular hydrogen. Results 
for the three sources of our sample are showed in Fig. 18. We 
found a rather constant value for this ratio, with weighted mean 
values of 0.94 + 0.11 for L1498, 1.14 + 0.11 for L1521E, and 
1.05+0. 10 for TMC 2, confirming the result achieved by Sarrasin 
et al. (2010) and Dumouchel et al. (2010), and consistent with 
detailed chemical models available in the literature (e.g. Herbst 
et al. 2000). 

This result resolves an important discrepancy between theory 
and observations which has lasted almost twenty years. There 
are rather few clear predictions of chemistry theory, but we can 
confirm that HNC and HCN seem to have similar abundances. 
However, given that what we measure is the column density of 
isotopologues, this implies that any fractionation should be the 
same for the two of them. 

Observations of higher-level transitions of HN'^C and 
H'^CN as wefl as H'^NC and HC'^N would help to refine the 
determination of column densities and would allow to estimate 
the isotopic [''*N]/['^N] ratio, too. 

8. Conclusions 

We have studied in this article the behaviour of the 7=1^0 
transitions of HCN(l-O), Hi3CN(l-0), and HNi3C(l-0) as a 
function of position in the three starless cores L1498, L1521E, 
and TMC 2. We also observed N2H+(l-0) and Ci^O(2-l) in 
TMC 2. Our main conclusions are as follows. 

1. Hi^CN(l-O) and HN'3C(l-0) are often assumed to be op- 
tically thin when computing column densities. Our results 
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Fig. 18. Abundance ratio [HNC]/[HCN] for the source sample as a function of the H2 column density. 



show that in the sources studied by us, this is inaccurate and 
indeed, optical depths are sufficiently high to make a reason- 
able estimate of the excitation temperatures of these transi- 
tions which can be compared with model predictions. 

2. The plot of H'^CN(l-O) excitation temperature against 
HN'''C(1-0) excitation temperature follows the curve ex- 
pected based on the collision rates recently computed by 
Sarrasin et al. (2010) and Dumouchel et al. (2010) thus con- 
firming these results. This plot also stresses the importance 
of calculations of potential surfaces and collisional coeffi- 
cients for isotopologues separately. Moreover these excita- 
tion temperatures coiTelate well with H2 column density es- 
timated on the basis of dust emission showing that dust emis- 
sion peaks really do trace peaks in the gas density. 

3. These latter results combined with our intensity-offset plots 
demonstrate convincingly that, at least in TMC 2, HCN and 
HNC survive in the gas phase at densities above 10^ cm""* 
where CO has depleted out. The implication of this is likely 
that CO survives at abundances of a few percent of its canon- 
ical value of around 10 "* and supplies the carbon required in 
lower abundance species. 

4. The profiles of the three satellites of HCN(l-O) become in- 
creasingly "skew" with increasing optical depth and indeed 
the corresponding H'''CN(l-0) profiles are reasonably sym- 
metric. This behaviour suggests the possibility of modelling 
the velocity field and abundance along the line of sight. 

5. We have used a model of the density distribution in L1498 
to describe the HN'^^C(l-O) and Hi3CN(l-0) results and 



find reasonable agreement with a model based on the pre- 
vious observations of Tafalla et al. (2006) containing a cen- 
tral "depletion hole" of radius 8 x 10'^ cm. This does not 
exclude models without a depletion hole, but does confirm 
our conclusions on the excitation discussed above. Indeed, 
rather surprisingly, our results suggest that HN'^C(l-O) and 
H'^CN(l-O) trace the high density nuclei of these cores 
when compared with other carbon-bearing species. 
6. Our results are consistent with the models of Herbst et 
al. (2000) who found that HNC and HCN should have simi- 
lar abundances in prestellar cores. 
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Appendix A: Non-LTE hyperf ine populations 

For a homogeneous slab with LTE between different hyperfine 
levels, one has 



1 - expi-fjT) 
1 - exp(-/;T) ' 



(A.l) 



where t is the total transition optical depth and /■ is the relative 
line strength of the /th component as in Table B.l. We applied 
this procedure to all the transitions, except for HN'^C(l-O), be- 
cause of the blending of the hyperfine components. Using the 
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same method adopted in Padovani et al. (2009) for C2H, we con- 
sider a homogeneous slab, then a two-layer model, and we com- 
pare different couples of ratios one versus the other to quantify 
the possible departure from LTE. As for C2H, we conclude that 
also a two-layer model cannot explain the observed intensities, 
requiring a proper non-LTE treatment. 

Comparing the upper and the lower panels of Fig. A.l re- 
lated to HCN(l-O) and H'^CN(l-O), respectively, it is clear that 
real departures from LTE are present and are stronger in L1498 
than in L1521E and TMC 2. In L1498, R32, the ratio between 
the weakest (88633 MHz) and the strongest (88631 MHz) com- 
ponent, is larger than unity and Rn > 1 in all the positions. 
This suggests a very high self (or foreground) absorption of the 
strongest component, as found for the C2H(l-0) emission in the 
same core (Padovani et al. 2009). The other two cores show ra- 
tios which deviate from the LTE curve revealing again the pres- 
ence of optical depth effects. As an instance, in the optically thin 
limit, i?32 should be equal to about 0.2, but the arithmetic mean 
of 7^32 is ~ 0.60 for L1521E and ~ 0.47 for TMC 2. 

Finally, an important remark follows from the H'^CN(l-O) 
ratios (lower panels): minor deviations from LTE are present 
even in this less abundant isotopologue which should be more 
optically thin. Even if optically thick lines are the foremost re- 
sponsible for the arising of non-LTE effects, the density can be 
high enough to allow transition within the same rotational level 
(e.g. the transition J,F - 1,2^ 1,1 has the same order of mag- 
nitude of rates relative to a transition between different rotational 
levels, see Monteiro & Stutzki, 1986). 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



Fig. A.l. Ratio of the integrated intensities of couples of components 
of HCN(l-O) (upper panels), and H''CN(l-0) (lower panels), where 
Rij represents the ratio between the integrated intensities W, and Wj 
(see Table B.l for component indices). Observational data: L1498 (yel- 
low squares), L1521E {green circles), and TMC 2 (black triangles). 
One-layer model (black solid curve), two-layer model (magenta dashed 
curves). 



Appendix B: Spectroscopic data and observational 
parameters 



Table B.l. HCN(l-O) and H"CN(l-0) frequencies of the hyper- 
fine components (From the JPL Molecular Spectroscopy Database: 
http : // spec . jpl . nasa . gov). 



Comp no. 


F' — F Frequency [MHz] 


/ 


HCN(l-O) 


1 


1-1 88630.4160 


0.333 


2 


2-1 88631.8470 


0.556 


3 


0-1 88633.9360 


0.111 


H'^^CN(l-O) 


1 


1-1 86338.7670 


0.333 


2 


2-1 86340.1840 


0.556 


3 


0-1 86342.2740 


0.111 



Table B.2. HN'^C(l-O) frequencies of the hyperfine components. For a 
complete list of the hyperfine components, see van Der Tak et al. (2009). 
In fact, there are four overlapped hyperfine components with = 2, 
three components with F'^ = 3, and three components with F = 1 — » 1. 



Comp no. 




transition 


Frequency [MHz] 


/ 






HN 


^'C(l-O) 




1 




- F, = - 


1 87090.675 


0.065 


2 




F^ = 2 


87090.791 


0.264 


3 




F^ = 3 


87090.834 


0.432 


4 




- F, = 1 - 


1 87090.886 


0.239 



Table B.3. Summary of observed molecules together with the observ- 
ing parameters: half power beamwidth, beam and forward efficiencies, 
system temperature, and precipitable water vapor. 



transition 


HPBW 


Bcff 


Feff 


T 

A sys 


pwv 




["] 






[K] 


[mm] 


HCN(l-O) 


28 


0.77 


0.95 


-130 


1-2 


H"CN(l-0) 


28 


0.77 


0.95 


-120 


1-2 


HN"C(l-0) 


28 


0.77 


0.95 


-140 


1-2 


N2H+(l-0) 


26 


0.77 


0.95 


-160 


1-2 


C'*'0(2-l) 


11 


0.55 


0.91 


-320 


1-2 
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